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O ■ ABSTRACT 

O ' 

We examine the feasibility of using dust-scattered X-rays for direct determination of 
distances to nearby galaxies with bright background AGNs, QSOs, or GRBs. We show 
how the Chandra X-Ray Observatory could be used to determine the distance to M31 
^ . to an unprecedented absolute accuracy of ~ 1%. 

Subject headings: scattering - X-rays: ISM - galaxies: individual (LMC, M31, M81) - 
£C) ' X-rays: general 
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1. Introduction 

o 



Or 



Distance measurement is one of the fundamental problems of astronomy. The distances to 
extragalactic objects are generally based on a "ladder": (1) Luminosities of nearby (typically low- 



2 ' luminosity) stars are calibrated using trigonometric parallax. (2) Luminosities of more luminous 

stars (e.g., Cepheids, or stars at the tip of the red giant branch) are then determined by observing 
such stars in apparent proximity to analogues of the trigonometrically-calibrated low luminosity 



stars. (3) Distances to nearby galaxies are then determined by photometry of what are believed 
to be analogues of the luminous stars. Errors are contributed by difficulties of photometric cali- 
bration, possible physical differences between the calibrated luminous stars and the extragalactic 
targets, possible blending, and errors in corrections for interstellar extinction. Bonanos et al. (2003) 
conclude that distances to galaxies in the Local Group are today known to no better than 10-15%. 

Geometric methods do exist for direct determination of extragalactic distances. Detached 
eclipsing binaries (DEBs) can be used (Paczyhski 1997), with distances to individual DEBs in the 
LMC and SMC claimed to be accurate to 4-14% (Guinan et al. 1998; Fitzpatrick et al. 2003; 
Harries et al. 2003). Bonanos et al. describe a project to determine the distance to M31 using 
DEBs. The "light echo" of fluorescent emission from the ring around SN 1987A has been used to 
estimate the distance to the LMC (Panagia et al. 1991; Gould &: Uza 1998 and references therein). 
Sparks (1997) has discussed the use of polarization in using light echoes to determine distances, and 
Xu et al. (1995) discuss using light scattered by dust to study the 3 dimensional structure of the 
ISM in front of SN 1987a. VLBI observations of the orbital motions of H2O masers in the galaxy 
NGC 4258 allow the distance to be determined to an accuracy of 4% (Herrnstein et al. 1999) if it 
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is assumed that the transverse motions of the maser spots coincide with the physical motion of the 
gas. 

Using time-delayed X-rays scattered by dust grains to estimate astronomical distances was 
originally proposed by Trumper & Schdnfelder (1973). Here we demonstrate the feasibility of 
determining extragalactic distances by repeated X-ray imaging of a time-varying AGN or QSO 
located behind a foreground galaxy. The distant point source will be surrounded by a time-varying 
X-ray halo. Measurement of the time-delay of the halo, as a function of halo angle, allows the 
distance to the foreground galaxy to be determined. Being based purely on geometry, the method 
is essentially free of systematic uncertainties. The question is whether the time delay can be 
accurately measured using a realistic observing program - this is the focus of the present study. 

The method depends on the X-ray scattering properties of interstellar dust grains, summarized 
in §2. The geometry of X-ray scattering by a foreground galaxy is described in §3. In §4 we consider 
some practical limits to the use of X-ray scattering to determine extragalactic limits. The method 
depends on variability of the X-ray source; in §5 we review the X-ray variability of AGNs, in §6 
their availability on the sky.. 

In §7 we present a general procedure for determining the distance to a galaxy from repeated 
images of a background AGN and a field around it. 

In §8 we discuss an observational program using the Chandra X-Ray Observatory to determine 
the distance to M31 using X-rays from the BL Lac object 5C 3.76. We carry out simulations 
to demonstrate the method, including photon-counting statistics and a realistic background. The 
reader interested primarily in the efficacy of this method may wish to go directly to Figure 8 and 
Table 1. A 5 Ms observational campaign with Chandra could determine the distance to M31 with 
an absolute uncertainty of ~4%, using only X-ray observations; a 10 Ms campaign can reduce the 
uncertainty to < 1%. We discuss how ancillary observations of CO and H I can be used to further 
improve the accuracy of the method. 

We consider using background AGNs to determine the distances to other galaxies in §9, but 
conclude that M31 is the most favorable opportunity. In §10 we assess the feasibility of using a 
background gamma-ray burst to determine extragalactic distances, but conclude that with present 
X-ray telescopes this method will likely be limited to the LMC and SMC. Use of X-ray telescopes 
other than Chandra is considered in §11. Our conclusions are summarized in §12. 

2. X-Ray Scattering by Dust 

Interstellar dust grains scatter X-rays through small angles, as was first pointed out by Over- 
beck (1965), Slysh (1969), and Hayakawa (1970). Because of this scattering, an image of an X-ray 
point source includes a "halo" of X-rays that have been scattered by dust grains along the line of 
sight. First observed by Catura (1983) using the Einstein observatory, scattered X-ray halos have 
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since been studied by a number of telescopes, including Einstein (Mauche & Gorenstein 1986), 
ROSAT (Predehl & Schmitt 1995) and Chandra (e.g., Smith, Edgar, & Shafer 2002). 

The scattered photons have a greater path length from source to observer than the unscattered 
photons. Therefore, if the source is time-variable, there will be a time-varying halo which is delayed 
relative to the observed variations in the point source. Because, for a given halo angle, the time delay 
depends on the distance to the source, Trumper & Schonfelder (1973) proposed that observations of 
time-variable scattered halos could be used to determine distances to variable X-ray sources. This 
effect has been used to estimate the distance to Cyg X-3 (Predehl, Burwitz, Paerels, & Trumper 
2000) and to Nova Cygni 1992 (Draine & Tan, in preparation). Very recently, Vaughan et al. (2004) 
used the X-ray halo around GRB 031203 - seen through dust with Ay rs 3 mag - to study the dust 
distribution toward t = 255 deg, b = — 4.6deg. 

There continue to be uncertainties about the composition, size, and geometry of interstellar 
dust grains, but observations of X-ray scattering by dust are in reasonable agreement with grain 
models that are approximately consistent with a broad range of observational constraints (Draine 
2003a). For a model of interstellar dust consisting of a size distribution of carbonaceous grains 
and silicate grains, the differential scattering cross section per H nucleon at X-ray energies can be 
approximated by (Draine 2003b) 



dn ) da 2iT9 s de s Li + (0 s /0 s ,5o) 

where cr sca is the total scattering cross section per H nucleon, 9 S is the scattering angle, and 

# s , 50 « 360" (^f) (2) 



V E 

is the median scattering angle for photons of energy E. For this model the dust has a total scattering 
cross section given by 

/ 771 \-l-8 

» 0.15 1- ) for 0.8keV < E < lOkeV . (3) 

Ay \ keV / 

where Ay is the V band extinction, in magnitudes. 



Scattering by Dust in a Foreground Galaxy 



Suppose that we observe a point source with flux F(t) located behind an intervening galaxy at 
distance D (see Fig. 1). Let the distance to the point source D s » D. Then X-ray photons observed 
at "halo angle" 9 have been scattered through an angle 9 S « 9/(1 — (3), where f3 = D/D s <^ 1. 

Consider a region in the galaxy subtending solid angle dQ, at angular separation 9 from the 
point source, with H column density N^. The differential flux contributed by this region is 



dH(t) 



(1-/3)2 



dcr sc . 
dQ* 



bg 



dn 



(4) 
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D s = D/jS (/5«1) 



B. T. Draine 3004-07.14.2301 

Fig. 1. — Geometry for single scattering of X-rays by dust in an intervening galaxy. 



where 5 is the increased light travel time for the scattered photons, and Ib g is background emission 
from unresolved sources and diffuse emission. For Euclidean space and 8 <C 1, the time delay for 
the scattered photons is 



s _ 1 D02 - 121 x 10?s ( D \ 
~ 1-/3 2c ~ 1-/3 V Mpc J V 100" 



(5) 



Now since 9 S = 9/(1 — (3), we can write 
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Consider annuli j, centered on the background point source, with inner radii subtending angles 



ipj, j = 1, 2, J. The flux from annulus j is 



Hj(t) = F(t-^)m° + I bg ,,% 



where 



(lpj+l/6h,50) 



is the "magnification" of annulus j, 



1 + (Vi+l/^,50) 2 1 + (^/^,50) 2 
''"scaj — -^Hj'Csca 



(8) 

0) 
(10) 



- 5 - 



is the X-ray scattering optical depth averaged over annulus j, 



is the r.m.s. value of 9 for the annulus, 

% = vr(^ 2 +1 - tf) (12) 

is the solid angle of the annulus, and ibgj is the background averaged over the annulus. The average 
time delay for annulus j is 

Note that portions of individual annuli can be "masked" - for example, if there are bright foreground 
point sources present in the image, or if a portion of the annulus falls beyond the detector boundary 
In such cases, the discussion below is unaffected provided only that the scattering optical depth 
T scai j and the background ibgj averaged over the annulus includes zero for the unusable "masked" 
regions of the annulus. It is straightforward to continue to allow the background to vary from one 
annulus to another, but henceforth we will assume the background Ib g to be independent of j. 

In the absence of other information, it is reasonable to assume the dust to resemble Milky Way 
dust (Draine 2003b), with scattering optical depth 

Tsca,j ~ A VJ x (^^j , (14) 

Here Ayj is the visual extinction averaged over annulus j, and T sca /Ay is given in eq. (3). 



4. Observability 

Let T bs be the time between the first and last observation. The largest useful halo angle is 

^(^^MW ' 

for halo angles 9 > max the time delay exceeds T ^ s . The radius R of the dusty region of the 
galaxy defines a second characteristic angle 

. (16) 



D V 5k P c / V D 

Let A6> be the angular resolution of the X-ray imager. Suppose that we wish to be able to determine 
the distance D to within a fractional error e, given a sufficient number of counts. There are 3 
different constraints which determine the maximum distance D for which this technique for distance 
determination is feasible. 
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i^yivi^v-.y B.T. Draine 2004.07.14.2301 

Fig. 2. — Characteristic angles as a function of distance D: Or is the angle subtended by 5 kpc. 6 mM (T b s ) is the 
angle corresponding to a time delay T b s . 650(E) is the median scattering angle for photon energy E. 



Even if the time delay Sj were measured precisely for some annulus j, accurate determination 
of the distance D = 2c5j/0'j requires measurement of the halo angle 0j to fractional accuracy 
e/2. Because the distribution of dust in the galaxy will be irregular, we cannot assume the 
dust to be uniformly distributed within an annulus. Thus determination of the scattering 
angle to a fractional accuracy e/2 requires an annulus width i^j+l ~ V'j < e %- The annulus 
width is at least equal to AO; thus AO < e0 max , or 



D < 



2e 2 cT obs 
(A0) 2 



0.827 Mpc ( 



V0.01 



(T ohs \ 
\W 7 s) 



A9J 



(17) 



2. The maximum halo angle is limited by the angular extent of the dusty portion of the galaxy. 
Measurement of halo angles to accuracy e would require AO < €0r, or 

eR / e \ / R \ ( 1" 



D < 



AO 



10 Mpc 



0.01 



5 kpc 



(18) 



3. Suppose r sca is uniform for r < R. The count rate of scattered photons then varies as 



[min(e max ,0ij)/0, 



h,50\ 



1 + [min(G max , R )/0 h:5O ] 



Ft 



(19) 
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Fig. 3. — Flux F relative to long-term average flux (F) for 5 independent realizations of light curves with variability 
amplitude A = 0.3 and correlation time r = 10 6 s, generated as described in Appendix A. 

From Figure 2 we see that for IMpc < D < 50Mpc and T obs < 10 7 s, 6 max <8r<. 6> Aj50 , so 
that the halo count rate H sca oc Ft SC3i /D, showing that the method is increasingly challenging 
as the distance D is increased. 



5. X-Ray Variability of AGNs 

We simulate the X-ray light curve of an AGN with the simple stochastic model described in 
Appendix A. The model light curve is controlled by three parameters: the average count rate (F), 
a characteristic correlation time r, and a variability amplitude parameter A. To choose realistic 
values of these parameters, we consider the X-ray lightcurves obtained for AGNs using EXOSAT 
(Lawrence & Papadakis 1993) and RXTE (Markowitz et al. 2003) and for QSOs using ROSAT 
(Ezoe et al. 2002). Based on these studies, we adopt a correlation time r « 10 6 s as representative 
(see Appendix A). 

The amplitude of variation can be characterized by F var , the square root of the variance 
normalized to the mean flux. Six Seyfert 1 galaxies observed by Markowitz et al. had values of -F va r 
ranging from 0.22 to 0.39, corresponding to the parameter A ranging from 0.22 to 0.38. We take 
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A = 0.3 as typical for AGNs. Five independent examples of synthetic light curves are shown in 
Figure 3. 



6. Density of Background AGNs 

The ROSAT North Ecliptic Pole Survey (Gioia et al. 2003), covering an 80.7 deg 2 area, detected 
54 AGNs with 0.5-2keV flux (corrected for absorption) greater than 2.5 x 10~ 13 erg cm -2 s" 1 . The 
catalog is approximately consistent with 

( dNiES ° >S) ) «1.6( inM S , X 312 ^ , (20) 

V dSl J E =ikeV V 10 13 ergcm 2 s 1 J 

where Se is the energy flux per unit photon energy E. Expressed in terms of count rate N for 
the back-illuminated CCDs of the Chandra Advanced CCD Imaging Spectrometer (ACIS-S BI) 
detector, we estimate 1 

m<N>N mm) / f, mm r 3/2 deg - 2 . (21) 

dn \ v .04cnts- 1 J 

If we have a foreground galaxy of area Af2, then we have a 50% probability of having at least one 
background source with a Chandra count rate exceeding 



i /l.6deg- 2 Af)\ 2/3 „ ,/A^ 2/3 

-r^- = .07cnts _1 5 

In 2 J Vdeg 2 



A> min = .04cnts- 1 ( • ^ o ) =.07cnts- 1 (- r - T ) . (22) 



7. Distance Determination 

Here we describe a simple method for estimating the distance from the observed X-ray images. 
Suppose that we have K images; image k covers "time bin" k extending from t± t k to i2,fc, with 

= t 2) k — ti,k- The images are time-ordered, with i 2 ,fc < ti,k+i- Each image is divided into 
circular annuli j = 1, J centered on the background source, with the innermost annulus exterior 
to the p.s.f. 

We now shift from discussion of fluxes to counts. Let Nj~ be the number of counts from the 
point source in time bin k. We assume that for each annulus j, the photon arrival times t arr are 
known; if not actually known, a random number generator is used to assign arrival times t arr for 
each detected photon within each integration interval. 



lr To estimate the Chandra ACIS-S BI count rates, we note that the OGLE source J005719.84- 722533.5 with 
(ES E )ikcV = 1.13 x 10 _13 ergcm _2 s _1 , seen through iV H = 1.5 x 10 22 cm~ 2 has a count rate (N) ~ 0.041 s _1 
(Dobrzycki et al. 2003). 
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The time delay for annulus j is Sj = a9?, with a an unknown parameter to be determined. 
There will be "background" counts due to noise in the detector, cosmic ray events, and unresolved 
celestial sources. Let I\, g A e gQj be the background count rate in the telescope in annulus j. We 
will suppose that this background I^ g (taken here to be uniform over the region of interest) is not 
known a priori, but needs to be estimated from the observations. 

For a trial value of the time delay parameter a, let the "exposure fraction" fj k (a) be the 
fraction of the time interval \t\^ + ctO?, t^^k + o/.O'j] during which the target was observed. Obviously 
fjk < 1 for t2 : k + o>9"j > t2,K, and fj k = for ti jk + a6j > ti,K\ fjk < 1 can also result from gaps 
in the observing campaign. 

Let Hjk(a) be the actual number of counts in annulus j during time interval [ti jk + ctQ?, t2,fc + 
a6j] (obviously, Hj k = when fj k = 0). From the Hj k , we construct the total (halo + background) 
counts in time bin k, 

J 

G k (a) = Y J H ]k {a) . (23) 
i=i 

For a trial value of the background ibg , we use the observations to estimate the "magnification" of 
each annulus, 

K 

m M,hg) = Wj(a)^2[H jk (a) - b jk (a,I hg )] , (24) 
k=l 

b jk (a,I hg ) = I hg A eS njf jk (a)At k , (25) 



Wj (a) 

For each trial (a, Ibg) we calculate 



K 



k=i 



(26) 



K 



v 2/_, r \-X^ w (n ) ~ M fe( Q ' V>jVfc - B k {a, I hg )} 
X (a,l hg ) = vV k {a) 27 — Tl ' ^ ' 

k= l a k \ a ^ I bg) 

where 

J 

B k (a,I hg ) = J2 b jk , (28) 

3=1 

J 

M k (a,I hg ) = ^ m j(^h g )f jk (a) , (29) 
i=i 

WM=°i&- . (30) 
22e=i G n«) 

The effective halo magnifications M k (a, i"b g ) depend on k because at early times the outer annuli 
are not usable (i.e., have fj k = 0) since the scattered light coming from them corresponds to the 
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unobserved point source light curve at t < t\^. The normalized weights W k , with ^ Wk = 1, are 
included to place greater weight on time bins where more counts have been observed. 

The denominator o\ in eq. (27) is calculated assuming photon counting statistics (see Appendix 

D): 

a 2 k (a, / bg ) = G fc (a) + N k M% - 2 (TV 2 + N 2 k M k + N k M k ) S k + (4iV 2 + N k ) S 2 k + 

(iV 2 + N k ) (T k + U k ) - 2N k V k , (31) 

where the estimated halo magnification M k (a, I bg ) is given by eq. (29), the weighting factors Wj(a) 
are given by eq. (26), and 

J 

S k (a,I hg ) = ^m jWj f 2 k , (32) 

3=1 

J K 

Tk(a,h g ) = ^2(wjfjk) 2 ^2{mjfjeN£ + bje) , (33) 

3=1 1=1 
K fj \ 2 

l=i \j=i J 
J 

Vk{a,I bg ) = ^wjfjkbjk ■ (35) 

3=1 

We require that fjk^k ^ 1 so that the statistical analysis of Appendix D be valid. 

There are K S> 1 halo epochs. We expect the function x 2 (o ; ,/bg) to have a minimum; let 
this be at (a*, /bg*)- If the parameters a and /bg are set to their true values, then (% 2 ) = 1 (see 
Appendix D). We adjust a and Ib g to minimize % 2 > but have a large number of annuli and epochs; 
we expect x 2 ( a *^bg*) ~ 1 if the only errors are due to Poisson statistics. Our best estimate for 
the distance is 

D* = 2(1 -(3)ca* . (36) 

While we do not claim the above methodology to be the optimal distance estimator - there may 
be alternative approaches that are less sensitive to Poisson noise - we demonstrate below by direct 
simulation that the method is capable of determining the distance for realistic data sets. 



8. Determination of the Distance to M31: Simulations 

8.1. Dusty Disk 

The dusty portion of the disk of M31 subtends a solid angle ~ 0.8 deg 2 (Xu h Helou 1996; 
Schmidtobreick et al. 2000), corresponding to a projected area = 1.10 x 10 45 cm 2 for an assumed 
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distance D « 770 kpc (Freedman & Madore 1990; van den Bergh 2000). The total dust and gas mass 
have been estimated to be (2.4 ±0.7) x 1O 7 M and M H « 2.5 x 10 9 M o within this disk (Xu & Helou 
1996, corrected to D = 770 kpc). Thus N u = 1.95 x 10 21 cur 2 . Using N u /A v = 1.87 x 10 21 cm' 2 , 
we estimate Ay ~ 1.0 averaged over the 0.8 deg 2 region. 

The disk thickness of ~ 200 pc is negligible compared to the distance D, so we can treat the 
disk as a scattering sheet. From eq. (5) we see that a time delay of 10 7 s corresponds to an angle 
9 w 100", so opposite sides of the scattering region will differ in distance to us by AD « 2D9smi ps 
0.001D sini, where i is the inclination. Inclination effects are therefore negligible, and the scattering 
region can be idealized as a thin sheet in the plane of the sky. 

8.2. 5C 3.76 

For fi = 0.8 deg 2 , eq. (21) would predict a 50% probability of at least one background source 
with Chandra count rate > 0.06 cnts -1 . Nature has provided at least 1 background AGN above 
this value: 2 5C 3.76 [= WSTB 37W051 = RX J0040. 2+4050; RA 00h40ml3.7s, DEC 40d50m05s 
(J2000)] - a compact nonthermal radio source about 40' (~ 9 kpc projected distance) from the 
center of M31. It was detected as an X-ray source by Einstein and ROSAT, and is classified as a 
BL Lac object (Perlman et al. 1996; Laurent-Muehleisen et al. 1999; Nilsson et al. 2003), of unknown 
redshift. 5C 3.76 is located behind iV(H I) = 1.34 x 10 21 cm~ 2 of atomic H (Dickey & Brinks 1993). 
If we assume that ~ 75% of the gas is atomic, with the balance either ionized or molecular, and 
assume that the dust is similar to Milky Way dust (with Nn/Ay ~ 1.87 x 10 21 cm -2 ), then we 
estimate Ay ~ 1.0 mag for this region of M31. 

Supper et al. (2001) reported a ROSAT PSPC 0.5-2 keV count rate of 0.124cnts _1 , corre- 
sponding to .ESV^lkeV) ss 3 x 10~ 12 erg cm -2 s _1 . 5C 3.76 was observed at 16 different epochs 
by the Chandra HRC, with a mean count rate 3 0.122 cnts' 1 (Williams et al. 2004). With 5C 3.76 
located behind iV(HI) ps 1.34 x 10 21 cm~ 2 , the X-ray spectrum is presumably strongly absorbed be- 
low 0.45 keV. If the flux were in 1.5 keV photons, for which the ACIS-S camera has ^cff ~ 600 cm 2 
(vs. ~ 200 cm 2 for the HRC) the ACIS-S on-axis count rate would be ~ O^cnts" 1 ; we adopt this 
as the estimated count rate. 

BL Lac-type objects are generally variable. The X-ray flux from 5C 3.76 was in fact observed 
to be variable by Williams et al.; the flux dropped by a factor ~ 1.7 between 2000-09-11 and 2000- 
10-12. Given its X-ray brightness, variability, and location, 5C 3.76 is an attractive background 



2 Somewhat outside the dusty disk, the background AGN Mrk 957 has a ROSAT count rate of 0.076 cnts -1 - 
about 60% of 5C 3.76 - and is known to show X-ray variability (Supper et al. 2001). Located just inside the D25 
countour of M31 (~20 kpc from the center), the dust extinction will be substantially lower than for 5C 3.76, making 
Mrk 957 less attractive than 5C 3.76 as a background source for determination of the distance to M31. 

3 5C 3.76 = object sl-75 was located 9.9 arcmin off-axis in the HRC observations. 
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source to use to determine the distance to M31. The variability spectrum of 5C 3.76 is not well- 
determined; we provisionally approximate it by our stochastic model with r = 10 6 s and A = 0.3 
(see Figure 3). 

Because the redshift of 5C 3.76 is not known, we do not know the value of (5 = D/D s . However, 
it is probably safe to provisionally assume a redshift z > 0.03, in which case D s > 125 Mpc, 
(3 < 0.006, and we can assume the factor (1 — /3) « 1 in eq. (36). Therefore even if the redshift is 
not known, the resulting uncertainty in the distance determination will be less than 0.6% provided 
z > 0.03. A program to determine the distance to M31 using X-rays from 5C 3.76 should, however, 
include optical spectroscopy of 5C 3.76 to attempt to determine its redshift. 

8.3. Diffuse Background 

Diffuse emission in M31 was measured by ROSAT (West et al. 1997). At ~ 40' from the 
center (along the major axis), unresolved sources contribute a ROSAT 0.5-2.0 keV count rate 
~ 7 x 10 -8 cnt s _1 arcsec -2 . Taking into account the greater effective area of the Chandra ACIS-S 
camera (600 cm 2 at 1.5 keV, vs ~ 150 cm 2 for ROSAT PSPC) we estimate a Chandra ACIS-S 
background count rate ~ 2.8 x 10 _7 cnts _1 arcsec~ 2 . 

XMM-Newton recently observed diffuse X-ray emission from a region 14' - 45' NE of the nucleus 
of M31, along the major axis (Trudolyubov et al. 2004), allowing an independent determination of 
the unresolved background. At ~ 30' from the nucleus, the diffuse 0.2-1.5 keV XMM count rate 
is ~ 9.2 x 10 _7 cnts~ 1 arcmin~ 2 (the lowest contour shown by Trudolyubov et al); from the radial 
variation seen by West et al. we estimate the count rate at 40' to be lower by a factor ~ 0.75. 
Taking into account the smaller effective area of Chandra 4 we estimate that the diffuse emission 
will contribute a Chandra count rate I® g A e f[ 3 x 10~ 7 cnt s _1 arcsec -2 , in excellent agreement 
with the background estimated from the ROSAT observations. 

8.4. Simulations 

With 5C 3.76 in mind, we simulate observations of a point source with mean 0.5 - 5 keV 
count rate (N) = 0.4 cnt s -1 , with a stochastic light curve generated following Appendix A with 
variability amplitude A = 0.3 and correlation time r = 10 6 s. We take the X-rays to have a 
characteristic energy E = 1.5 keV, and approximate the disk of M31 as a uniform scattering screen 
with a scattering optical depth r sca = 0.072 for the X-ray photons (corresponding to Ay = 1 mag). 
The scattering phase function for the X-rays is assumed to be that estimated for Milky Way dust 



4 At 1.5 keV, Chandra ACIS-S and XMM have effective areas of ~ 600 and ~ 1300 cm 2 , respectively. 
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Fig. 4. — Symbols: estimated values of rrij vs. annulus number j for a randomly-selected simulation with T b s = 
7 x 10 6 sec, for the best-fit distance D and background Ib g - Solid line: m° used in the simulation, rrij is largest for 
large j because the outer annuli (with AO = 1") have larger area. 

at 1.5 keV (see eq. 1,2). 5 The time delay parameter is set to Qo = Dq/(2c), where Dq = 770 kpc 
[and we take (3 = in eq. (13)]. 

In addition to the scattered X-rays, we assume a uniform background with a count rate 
^bg^cff = 3 x 1CT 7 cnt s _1 arcsec~ 2 . 

For simplicity, we simulate observing campaigns without gaps, with continuous exposures 
?obs = 3, 5, 7, and 10 Ms. For each case we perform 10 3 independent simulations. 

For the 10 Ms observing campaign, the maximum useful scattering angle [see eq. (15)] is G max = 
104". For 1.5 keV photons, the median scattering angle is #50 = 240". The mean count rate for 
scattered photons within 100" of the point source will be ~ 0.4-0.072-(100/240) 2 /[l + (100/240) 2 ] = 
0.0043 cnt s _1 , only ~45% of the background count rate within this region. 



5 While the total scattering cross section [eq. (3)] varies as ~ B -1-8 , the differential scattering cross section at 
< 83,50 [see eq. (2)] is nearly independent of E (see Fig. 8 of Draine 2003b). The observations considered here are 
at 9 < # s ,5o, and thus our estimated count rates for scattered photons are insensitive to the assumed energy spectrum 
of the point source. 
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Fig. 5. — Effective total halo magnification Mk for time bin k. Solid line: M° for the simulations. Symbols: Mk 
for three randomly-selected simulations with T b s = 7 x f 6 s. 

We use uniform time bins of width At = 5 x 10 4 s, so that the point-source counts per time 
bin (N%) « 2 X 10 4 . The halo annuli are constructed as described in Appendix B, with minimum 
annular width of A9 = 1", giving J = 101 annuli within 104". 

For each simulated data set, we construct the delay-corrected halo light curve Gk(ct) for a trial 
value of the time delay parameter a. For different trial values of I\j S , we calculate x 2 ( a i -^bg) using 
eq. (27). We repeat the procedure for different trial values of a, and adopt the values (a*, ibg*) f° r 
which x 2 i s minimized. 

For small values of T Q b s , the number of halo counts is small, Poisson fluctuations are very 
substantial, and a+ may be far from the "true" value «o- Because the amount of useful data 
decreases as the trial value of a (i.e., D) is increased, Poisson noise may result in a false minimum 
of x 2 for large values of D. Because we have a-priori limits on the plausible range of distances to 
M31, we search for the minimum of x 2 only over the range [0.78-Do, 1.22Z?o]- In the event that this 
minimum occurs outside the range [0.8Dq, 1.2Dq], we reject the distance determination altogether. 
Such rejections are, however, relatively rare, especially for the longer simulated campaigns. In a 
real observing campaign, the observations would be continued to accumulate data and eliminate 
such spurious results. For 5 Ms exposures, we had only 27 rejections in 1000 simulations; for the 
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Fig. 6. — (Gk — Bk) (black) and NkMk (red) for four different trial values of a, for a simulation with T b s = 5Ms. 
At late times (Gk — Bk) and NkMk go to zero because only the innermost annuli provide data. To the eye, NkMk 
and (Gk — Bk) appear to be more similar for the true delay a = ao than for the trial values 0.85ao, l.lao, or 1.15ao, 
but the Poisson noise is large. 

10 Ms campaign, we had rejections in 1000 simulations. 

The method involves estimation of the annular magnifications rrij from the observations. Figure 
4 shows the rrij estimated from a 7 Ms simulation. Because of Poisson statistics, the rrij are quite 
noisy. However, the (see eq. 29 and Figure 5) are much better behaved, as expected. 

Our method seeks the time delay that minimizes the difference between (Gk — Bk) and NkMk- 
Figure 6 shows these two functions of time for one T Q b s = 5Ms simulation, for four trial values of a. 
To the eye, N^M^ and (Gk — Bk) do appear to be more similar for the true delay a = a® than for 
the trial values 0.85ao, l.lao, or 1.15ao, but the Poisson noise is obviously large. Our x 2 statistic 
is intended to quantify the difference between (Gk — Bk) and NkMk, where (x 2 ) = 1 for the true 
time delay. 

In Figure 7 we show x 2 vs D/Dq for different trial values of I^ g , for 4 independent simulations 
with T Q bs = 7 x 10 6 s. The value of ibg giving the x 2 (a) with the lowest minimum ranges from 
0.87/^ to 1.092]^ - the background estimate is accurate to about 10%. It is also apparent that 
the location of the minimum is not highly sensitive to the estimate of 2b g . Most importantly, it is 



-16- 




Fig. 7. — x 2 vs D/Do for 4 randomly selected simulations with T b s — 7 x 10 6 s. In each case we show X 2 ( Q ) f° r 3 
trial values of Jb g : 0.927° g , 1.087° g , and the background that leads to the lowest minimum of X 2 { a ) F° r each case we 
give the estimate -Zbg,* as a fraction of I° g . 



seen that the minimum a* is very close to the true value ao - within ~3% - for each of these four 
simulations. 

Histograms of D^/Dq are shown in Figure 8 for each of the 4 values of T b s - Table 1 gives 
the median, 68% confidence interval, and rejection fraction, for 4 different values of T ^ s . For the 
shortest observing campaign, T Q ^ S = 3 Ms, the simulations result in a fairly broad distribution 
of derived distances D, with about 83 of the experiments resulting in "rejection". For the 917 
nonrejected cases, the median D/Dq = 1.026, and the 68% confidence interval extends from about 
0.967 < D/Dq < 1.137, i.e., "one cr" uncertainties of t^S"- 

As the observing time T G b s increases, the amount of useful data increases as ~ T^ hs , because 
it is possible to use scattered halo photons at larger angular separations from the point source. As 
a result, the accuracy of the the distance determination improves rapidly. As seen from Figure 8 
and Table 1, it appears to be possible to determine the absolute distance to M31 to within ±1% 
with a 10 Ms observing campaign with Chandra. 
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Fig. 8. — Histograms showing distribution of D CBt /Do for M31, where D es t is the estimated distance, for simulations 
with T bs = 3, 5, 7, and 10 Ms. In each case, we show results for 10 3 simulations. Distance estimates differing from 
the true distance by more than 20% are rejected (see text). Of the remaining distance estimates, 68% fall between 
the dotted lines, 16% above, and 16% below; the dotted lines can be loosely understood as ±la uncertainties in the 
distance estimation. 



8.5. Discontinuous Observing 

Our discussion assumed, for simplicity, continuous observations of M31. Realistic orbital and 
scheduling considerations would require gaps in the observing. The methodology presented in §7 
can be applied to observing campaigns that include gaps. 
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Table 1: Results for Simulated M31 Observing Campaigns 



Tobs 


D cst /D a 


rejects 


3 Ms 


1-0261;^ 


8.3% 


5 Ms 


1.007^ 


2.7% 


7 Ms 


1.005l;gg 


2.1% 


10 Ms 


l-002±;ggg 


< 0.1% 



a 68% confidence interval based on 10 3 independent simulations (see text), rejecting simulations where the distance 
estimate differs from the true distance by more than 20%. 

The simplest approach would be to use only the time intervals for which the AGN was observed, 
and only the counts from annuli j with arrival times i arr such that t arr — a9j was a time when the 
AGN was observed. This approach would, however, have the disadvantage of discarding halo counts 
for which i arr — a6j falls in a gap when the AGN was not observed. 

A better approach, if individual gaps are not too long, would be to interpolate to estimate the 
AGN counts Nk for time intervals falling in the gaps in the observations. This can be accomplished 
using the procedure described by Press, Rybicki, h Hewitt (1992). Using these interpolated counts 
in the same way as the actual counts Nk, one could use the same procedure described above to 
minimize x 2 { a ,Ibg) to find the distance D. 

Provided the interruptions in observing are sufficiently brief (< 1 day) that the AGN lightcurve 
can be reliably interpolated, discontinuous observing should be almost as effective as continuous 
observations with the same total exposure time. 

8.6. Using Multiwavelength Data 

We have shown above that it is possible to determine the distance to M31 using only X-ray 
imaging. In particular, the X-ray data alone is used to estimate the "magnification" rrij of each 
annulus, which is proportional to the amount of interstellar dust present in the annulus. 

Imaging at other wavelengths can also provide information on the distribution of interstellar 
dust. Note that the scattering zones j need not be complete annuli - the annuli can be subdivided 
in the event that the dust distribution is nonuniform. Indeed, if a portion of the annulus were 
known to be dust-free, this should be masked off, as it would contribute Poisson noise from the 
background but no signal. If the dust distribution over the annulus is appreciably nonuniform, it 
may be advantageous to subdivide the annulus into subannuli, with low and high dust densities. If 
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H I and CO observations are used to estimate the surface density of gas, the magnification 

rrij = f(9) (jlp^J % [X H X,(H21cm) + 2X C 0^(CO 2.6mm)] (37) 

where f(9) is the angular distribution function for X-ray scattering 6 at halo angle 9, Qj is the solid 
angle of (sub)annulus j, ij(line) is the line intensity in a spectral line averaged over (sub)annulus 
j, Xh is the conversion factor from 21cm line intensity to H I column density, and Xqo is the 
conversion factor from CO line intensity to H2 column density. 

There are at least two different ways in which this information could be used: 

1. If the angular resolution and signal/noise ratio of the 21cm and CO images is sufficient, and a 
good estimate for f(9) is in hand, we could use eq.(37) [rather than eq. (24)] to estimate the 
rrij for the (sub)annuli, but with the scaling factor t scs ,/Nu to be determined. Now instead 
of trial values of a and I^ g , we would instead consider trial values of a and T sca /Nn; for each 
(a, T sca /./VH) pair, the background Ibg is estimated from 

4 g (a, r sca /iV H ) = V L (38) 

A cS Z^fc=l iV fc 1^3=1 U JJjk 

Using this value of I\, g , x 2 («> ^sea/Mi) is then calculated for the trial values of a and r sca /A^H- 
The best estimates of a and T sca /iVH are then found by minimizing x 2 ( a > ^sea/Mi)- 

2. Alternately, the magnifications rrij(a,I\ )g ) estimated from the X-ray data alone [eq. (24)] 
could be compared to the rrij estimated from eq. (37), by computing an error function 

2 = [mj (eq.24) - (r scll /N H )sj)} 2 ^ 
' a 2 (mj) 



j 



Sj = f(9)£lj [X H ij(H21cm) + 2X co Ij(CO 2.6mm)] (40) 



Tsca Ef=i m j("> Ib g )sj/(J 2 {m^ 



■3) 



(41) 



^ e;=i^ 2 /^k) 

a 2 (rrij) is the variance in rrij due to Poisson statistics [see eq. (D23)]: 

a 2 (m,) « rrij ( Wj f £ f jk N k + (^-) 2 £ 6 jfc + (nijiu,) 2 £ /? fc JV* . (42) 

k k k 

For the wrong estimate of the distance, the estimated magnifications rrij will be systematically 
in error. Thus Xd( a ) should have a minimum at the true value of the time delay parameter 

Q 



6 With normalization / f(6)2n sin 6d6 = (1 - /?)" 
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For method (1) to be useful, it is necessary to have imaging data at angular resolutions smaller 
than the larger of (a) the annular widths (1")) or (b) the angular scale over which the dust surface 
density varies. It would be useful to have 21cm and CO interferometric imaging of this 4' x 4' region 
of M31. In principle, the dust could be observed directly using SCUBA on the JCMT (HPBW 
of 7.5" at 450/xm), or MIPS on the Spitzer Space Telescope (~ 30" at 160/xm), but the angular 
resolution of these instruments is probably insufficient to trace the dust variation from annulus to 
annulus. 

The best use of multiwavelength data will depend on its angular resolution and signal/noise 
ratio. We do not undertake any such simulations here. However, it is clear that this additional 
information can only improve the distance determinations. If an observing campaign is undertaken 
with Chandra, use of multiwavelength data should be further investigated to make best use of the 
Chandra observing time. 

9. Other Galaxies: LMC, SMC, and M81 

With an H I mass 3.1 ±0.6 x 1O 8 M (Luks & Rohlfs 1992), the LMC has (JV H > = 1 x 10 21 cm" 2 
over f2 = 40deg 2 . The dust/gas ratio in the LMC is ~ 40% of the Milky Way value, hence 
we estimate {Ay) k, 0.21 mag, and r sca « 0.015(^/1.5 keV)" LS . What is required is a bright 
background source. Dobrzycki et al. (2003) report 5 new X-ray QSOs behind the SMC, the brightest 
of which has a Chandra count rate 0.041 cnts -1 . Over the 40deg 2 area of the LMC, eq. (22) leads 
us to expect the brightest background AGN to have ESe ~ 2x 10 -12 erg cm -2 s -1 , with a Chandra 
count rate N ps 0.8 cnts -1 . This is only twice the count rate for 5C 3.79, and r sca is likely to 
be a factor of two lower than for 5C 3.79, so the actual halo brightness for an AGN behind the 
LMC may not be higher than for 5C 3.79 behind M31. The reduced distance (D = 66kpc for the 
LMC) allows observation of the time-delayed halo out to ~ 112"(T o b s /10 6 sec) 2 , but an accurate 
distance determination would still require a substantial observing campaign, and ultimately would 
be limited by the uncertain three-dimensional geometry of the LMC. In many ways M31 is a more 
attractive candidate for X-ray distance determination using AGNs. 

As discussed in §4, distance determination using X-ray halos can be applied to galaxies beyond 
M31, but with increasing difficulty. M81, at a distance D = 1.4 Mpc, is perhaps the next most 
likely candidate. Immler & Wang (2001) present ROSAT observations of M81. Aside from the 
nucleus of M31, the brightest source in the field is source X9, located 12.2 arcmin E of the nucleus 
(5 kpc projected separation). Variable (with an amplitude exceeding a factor of 2.5; Immler & 
Wang 2001), with a ROSAT PSPC count rate of 0.20 cnts" 1 , and located in a region of M81 with 
iV(H I) ~ 3 x 10 21 cm" 2 (Immler &: Wang 2001), X9 would be suitable for distance determination 
if it were a distant background source, but Wang (2002) argues that X9 is an intermediate-mass 
black hole associated with M81. At this time there are no known bright background sources that 
would allow X-ray distance determination to M81 with existing observational capabilities. 
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10. Using X-Rays from Gamma-Ray Bursts 

Suppose a gamma-ray burst (GRB) occurs behind a nearby galaxy (e.g., LMC, SMC, or M31), 
with fluence Fq in 1-2 keV X-rays. We approximate the GRB X-rays as a single short-duration 
pulse. The galaxy is approximated as a uniform thin screen with X-ray scattering optical depth 
r sca . The dust is assumed to have the scattering phase function of Milky Way dust for 1.5 keV 
X-rays [see eq. (1,2)]. 

The fluence in scattered X-rays is 

dF SCSL = F T sca —^2iTed6 (43) 



? 

■[l + (0/0 h ,5Q 

Since the geometric time delay t = a9 2 , we find 



^0Tsca ri d( ^ 5 °\ 21 2 / ( 44 ) 



dF sc& = F r sca ' , (45) 

[1 + r/th,5oJ 

where 

*h,50 = a9l >50 = 3.9 x 10 6 s (^^) • (46) 

We now suppose that we image over the interval (t — At/2, t + At/2). The scattered photons will 
be observed as a rin g with radius 9 = {t/a) 1/2 and width (A6% bs = max[A0, At/(2a9)], where A9 
is the angular resolution of the X-ray imager. The number of scattered photons counted in the ring 
will be 

N r « F WU f ^ (47) 

(1 + V*h,50) 

„ / Fp \ / r sca \ /56kpc\ /240^\ 2 / A cS \ At/10 5 s 

~ Vl0 2 cm-2; 1 0.015 J\ D J\9 h:50 ) \m0cm 2 J (1 + t/t h ) 2 ' 1 ° j 

where A e s is the effective area of the X-ray imager. 

The ring will therefore be detected with a signal-to-noise ratio 

SIN « — , 1/0 , (49) 

1 (iV r + A r bg )i/2 ^ > 

where the number of background counts in the ring is 

N hg = I hg A eS 2Tr9(A9) ohs At . (50) 

As an example, consider a GRB with a 0.5 — 5keV X-ray fluence Fq = 10 3 photon cm~ 2 , 
located behind a region of the LMC (D = 56kpc) with Ay = 0.2 (r sca « 0.015). A At = 2 x 10 4 s 
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exposure with the Chandra ACIS (A = 600 cm 2 ), taken at a time t = 10 6 s (when the ring radius 
is = 122") will have AO = 1.22", and will detect N r = 29 counts due to the scattered halo. If 
-^bg^efr = 3 x 10~ 7 cnt s _1 arcsec -2 , there will be iVbg ~ 6 background counts in the annulus, and 
the ring will be detected with a signal-to-noise ratio S/N ~ 5. With a width AO = 1.22", the 
average ring radius, and therefore the distance D, could be determined to better than 1%. 

The above example assumed a relatively high X-ray fluence Fq = 10 3 photon cm -2 ; if the 
fluence is lowered to Fo = 10 2 photon cm -2 , the ring would have a signal-to-noise ratio of only 
0.5 in a single At = 2 x 10 4 s exposure. However, by taking, say, 50 such exposures in succession 
and appropriately stacking them to allow for expansion of the ring from one exposure to the next, 
one could increase the S/N ratio to 0.5^/50 = 3.5. RXTE, with a FOV of 1130 deg 2 , detected 
8 GRBs in the last 9 months of 1996 (Smith et al. 2002) corresponding to an all-sky event rate 
of ~ 320yr _1 . The detected sources had a median 1.5-12 keV fluence of ~ 10 2 photon cm~ 2 , so 
we estimate an all-sky event rate of ~ 10 2 yr _1 of GRBs with Fb(0.5 — 5keV) > 10 2 photon cm~ 2 . 
Since the LMC has a projected area of 40 deg 2 , the probability per unit time of a GRB with 
F (0.5 - 5keV) > 10 2 photon cm" 2 behind the LMC is ~ O.lyr" 1 . 

While GRBs could also be used for distance determination for M31 and other spirals, the 
ring counts N r oc 1/D (see eq. (48), and the much smaller angular size (0.8 deg 2 for M31, vs. 
40 deg 2 for the LMC) makes it very unlikely to have a suitably bright GRB behind M31. Distance 
determination using X-ray halos around GRBs will likely be limited to the LMC and SMC. 

11. Other Telescopes 

XMM-Newton offers a larger aperture than Chandra, but its relatively poor angular resolution 
(AO ~ 6") makes it unsuitable for precision distance determination using scattered X-rays. The 
proposed Constellation X mission is projected to provide an order-of-magnitude increase in collect- 
ing area (A e g « 13000cm 2 @ 1.5 keV), but current plans call for angular resolution AO « 6", which 
would probably preclude accurate distance determination. If the angular resolution were improved 
to ~ l", Constellation X would be capable of extragalactic distance determination to local group 
galaxies using background AGNs with ESe ~ 10~ 13 ergcm~ 2 s -1 . 

12. Discussion and Summary 

Current distance uncertainties to M31 are of order 10%. Bonanos et al. (2003) are undertaking 
to use detached eclipsing binaries in M31 to establish the distance with an anticipated accuracy of 
5%. We have shown above that a Chandra program of extended observations of the field around the 
background X-ray source 5C 3.76 has the potential to allow the distance to M31 to be determined 
to an absolute accuracy of ~ 1%, for a 10 Ms observing campaign. In the course of an observing 
campaign, we estimate a 90% probability that a distance estimate accurate to ~ ±10% would 
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emerge after 3 Ms (see Table 1); the distance estimate continues to improve as more data accumu- 
lates. We caution that our observing time estimates are based on assumed variability properties 
for 5C 3.76; further observations to characterize its variability should be undertaken. In addition, 
monitoring of 5C 3.76 would allow a distance-determination campaign to be "triggered" to begin 
at a time when the source undergoes an outburst (as opposed to our simulations, which started at 
random times). 

A 10 Ms observing campaign would constitute a major commitment of Chandra observing 
time, but would allow determination of the distance to M31 to unprecendented absolute accuracy. 
There do not appear to be any other methods capable of determining extragalactic distances to 
such precision in the foreseeable future. A precision determination of the distance to M31 would be 
of great importance to calibrate the luminosities of stars, and to allow accurate determination of 
the Hubble constant. The proposed observing campaign would, of course, also provide a census of 
the X-ray source population in this region of M31 down to lower luminosities than current surveys, 
as well as information on variability of the point sources. 

To summarize our principal conclusions: 

1. Realistic X-ray scattering properties for Milky Way interstellar dust have been used to esti- 
mate the strength of X-ray scattering halos around distant X-ray point sources seen through 
foreground galaxies. 

2. The X-ray variability properties of AGNs and QSOs are reviewed, and a method for generating 
simulated light curves is developed (Appendix A). 

3. We present a simple method to use X-ray observations of a background AGN and the scat- 
tered X-ray halo around it to directly determine the distance to a foreground galaxy. The 
methodology allows discontinuous observations to be employed. Provided that the fraction 
of the time lost to the gaps is not large, discontinuous observations will be almost as effective 
as the idealized continuous observations simulated here. 

4. We use simulations to demonstrate the feasibility of direct determination of the distance to 
M31 using the background BL Lac object 5C 3.76. For the variability properties that we 
have assumed, a 5 Ms observing campaign with Chandra would allow the distance to M31 to 
be determined to ~ 4% accuracy, and the absolute distance uncertainty could be reduced to 
<1% with a 10 Ms campaign (see Figure 8 and Table 1). 

5. We discuss how observations at other wavelengths -HI 21cm and CO 2.6mm aperture 
synthesis imaging, in particular - can be used to complement the Chandra observations. 

6. We consider distance determination using the X-ray halo around GRBs. This is a viable 
method for determining the distance to the LMC or SMC, although GRB statistics suggest 
that sufficiently bright GRBs occur behind the LMC at a rate of only ~ O.lyr -1 . 
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A. Synthetic AGN Light Curve 

We describe a simple procedure to generate a stochastic light curve with statistical properties 
resembling X-ray light curves for AGNs. 

Let gj be a sequence of independent random variables with zero mean and unit variance: 
(<7i) = and {gigj) = For times tj = to + jh, let 

v 3 = (1 - e" 2 Y /2 £ me-*-* ■ (Al) 

i=— oo 

If we set A = h/r, then 

( Vi ) = , (A2) 

{viVj ) = e -\U-t 3 \/r . (A3) 

Consider now a (positive-definite) flux given by 

F(tj) = ae Av * (a > 0) . (A4) 

If the gj are assumed to be gaussian random variables, then 

(gj) = for n odd, (A5) 
n\ 

= — T7T, — ttt for n even; (A6) 
2"/2( n /2)! v ' 

for this case one can show that 

(F n ) = a n e nlA2 l 2 for n = 0,1,2,3,... (A7) 
The r.m.s. fractional variance 



F = 

1 var — 



[<F 2 > - (F) 



21 1/2 



so that the coefficient A can be obtained from F var 

?2 Ml/2 



A 2 

e A - 1 



1/2 



(A8) 



A=[ln(l + F 1 2 ar )] i ^ . (A9) 
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Fig. 9. — Structure function and PSD for synthetic light curves for different values of A. 



From (A7) and (A9) we see that the coefficient a in (A4) is 



a = (F)e" A2 / 2 



(F) 



(1 + ^ar) 1/2 



The variability can be described by the structure function 

SF(r)^([F(t)-F(t-r)] 2 ) 

or the power spectral density 

PSD(/) ee ^|F(/)| 2 , 

J- obs 



(A10) 



(All) 



(A12) 



where F(f) is the discrete fourier transform of F(t), T Q b s is the total time span observed, and At 
is the time between measurements. 

In Figure 9 we show the structure function and power spectral density obtained from these 
light curves. The dimensionless quantity / x PSD(/)/(F) 2 peaks at f p ps 1/2ttt. For A = 0.3 the 
maximum of / x PSD(/)/(i ? ) 2 is ~ 0.015; this is in agreement with the maxima of / x PSD(/) for 
the 6 Seyfert galaxies observed by Markowitz et al. The peak frequency f p < 10~ 7 Hz for Fairall 9, 
NGC 5547 and NGC 3516, ~ 10~ 7 Hz for NGC 4151, and ~ 10~ 6 Hz for Ark 564 and NGC 3783. 



B. Choice of Annuli 



Ideally, annuli would be sufficiently narrow that the variation in time delay across an annulus 
would be a fraction e <C 1 of the time r characterizing fluctuations of the X-ray source. However, 
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the minimum annulus thickness is set by the angular resolution A6 of the X-ray imager. Thus we 
take the outer boundaries to be 



i>j = V>iJ 1/2 for j < j c 

= te 1/2 + U ~ 3c) AO for j>j c 



Jc 



nint 



ecr 



nint 



2.07 



A9 



[2D(A8Y\ 

where nint(x) is the integer nearest to x. 

Note that annuli j < j c are uniformly spaced in time delay 8. 
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C. Simulation 

Our simulated observating campaigns were created as follows. We adopt a value for the time- 
averaged AGN count rate (F)A e g, the AGN correlation time r, the total exposure time T Q b s , the 
scattering optical depth r (assumed uniform), and the time delay coefficient a, where the extra 
light travel time for a scattered photon from angle 6 is a6 2 . 

We assume the observing campaign to capture K images containing the AGN, each for an 
exposure time T ^, s /K, where K = 100. For simplicity, we assume that there are no time gaps 
between the images. For an assumed value of the time-averaged point source count rate (F) and 
correlation time r, we compute a random realization of the AGN light curve F(t) using eq. (A4), 
using short time steps h Ct. 

The simulated images are produced as follows: The point source counts are calculated by 
randomly drawing from a Poisson distribution with expectation value J F(t)dt integrated over the 
exposure k. 

Because the time delay can vary significantly across a single annulus, we divide each annulus j 
into a large number of subannuli. Each subannulus has a magnification calculated using equation 
(9), with the boundaries ipj replaced by the boundaries of the subannulus. We assume the back- 
ground ibg to be uniform across the image. Each subannulus has an expected number of counts 
during the exposure time Ai&; we sum to obtain Hj k , the expected number of counts for annulus j 
during exposure k. We then draw the "observed" counts Hjk randomly from a Poisson distribution 
with expectation value Hj k . 

The simulated counts and Hj/, now contain Poisson "noise" . We take and Hj^ and use 
the method described in §7 to estimate the background Ibg and time delay coefficient a. 
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D. Statistical Errors 

Let raj be the true magnification of annulus j, let the point source have expected number of 
counts N® in time bin k, and let the true background intensity be . Recall that the exposure 
fraction fj k is the fraction of the time At k for which annulus j was observed. 

Suppose we have observed counts Hj k in annulus j and time bin k. For the true time delay 
parameter a , we proceed (as discussed in §7) to construct the halo light curve G k (ao). Let 

A k (a) = A cS Y,^]fjk{a)At k . (Dl) 

j 

For each time bin k, the observed G k (ao) differs from the expected value [M®N% + 1® A}.(ao)] due 
to Poisson fluctuations. The expectation values iV° are not known, but the N k are observed. The 
true annular magnifications are not known; the estimated annular magnifications mj(ao, I^ g ) 
and cumulative magnifications M k (ao, ) are obtained from the observed N k and the observed 
G fc (ao) using eq. (24, 29). 

Assuming that we have guessed the correct value of Ibg = ^b g > we wish to calculate the 
expectation value ([G k (a ) - lg g A k (a ) - M k (a , lg g )N k } 2 } iorl<k<K. Let 

J 

M k = £™?/i*(ao) (D2) 

b% = I° g A eS njf jk (a )At k (D3) 

Then the observed N k and Hj k are 

N k = N° k +u k (D4) 
H jk = m° j f jk N% + b% + x jk (D5) 

where v k and Xj k are independent random variables. The N k and Hj k obey Poisson statistics, thus 

H> = (x 3k )=0 (D6) 

<^ 2 > = {ul) = Nl (D7) 

<^ 4 > = 3K) 2 + iV fc (D8) 

(x%) = m]f jk Nl + b° jk (D9) 

The magnifications rrij (ao, I® g ) are estimated using eq. (24). Expanding in powers of EfcLi fjkN®] -1 , 

sr^K tt h o K K 



7 a/" w m ? + w ? E - m >^° E ^ ' (D10) 



"'j i- "j 

k=l k=l 

where 

„o 



K 



k=i 



ocl . (Dll) 
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We find 



where 



([G k -I° g A k -M k N k Y) « N° k M k (l + M k ) + B° k 

-2 [(N° k ) 2 + {N Q k fMl + iV°M fc ] 5 fc ° 
+ [4(iV fc °) 2 + Nl] {Si? 

+ [(N° k ) 2 + N° k ] (T° k + U° k ) - 2N° k V k ° , (D12) 



B° k = E4 (° 13 ) 

J k 



i! - E «fe) 2 E Kwv? + «&) ~ ^ + f) ~ w + w (D15 > 

3=1 l=\ \ / 

U° k - E^(E-K/^) ~iv(^)^M (D16) 



l=i \j=i 
J 



V* - E w °J^ b %-YN (D17) 

where we have indicated the dependences on M k , N = J2 k N k , and B = Y^ fc B k . Eq. (D12) is based 
on an expansion in powers of (Yl k fjk^k) 1 ~ J/N, and therefore is applicable only if individual 
annuli are large enough that J2k fjkN k S> 1. Since M k <C 1, in eq. (D12), the terms in (S k ) 2 and 
U k are normally negligible. The leading order terms are 

KM* ~ ^ (D18) 



B° ~ | (D19) 



(NlfSl ~ ^ (D20) 

(*2W ~ ^ + ^ (D21) 

A« ~ A (D22) 

It is also useful to estimate the dispersion of the estimated magnifications rrij around the true 
values m°y From eq. (D10): 

((m 3 - m ») 2 ) « m; K) 2 E + K) 2 y: b% + KS°) 2 E 4*2 • (° 23 ) 

A; A; fe 
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